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ABSTRACT 


The  frequency  (or  phase  velocity)  of  axially  symmetric  free 
vibrations  in  an  elastic,  isotropic  circular  cylinder  of  medium 
thickness  is  studied  on  the  basis  of  the  three-dimensional  linear 
theory  of  elasticity  and  also  on  the  basis  of  several  different 
shell  theories.  To  be  in  good  agreement  with  the  solution  of  the 
three-dimensional  equations  for  short  wave  lengths,  an  approximate 
theory  has  to  include  the  Influence  of  rotatory  inertia  and  transverse 
shear  deformation  for  example  in  a  manner  similar  as  in  Mindlin’s 
plate  theory,  A  shell  theory  of  this  (Timoshenko)  type  is  deduced 
from  the  three-dimensional  elasticity  theory.  From  a  comparison  of 
phase  velocities  it  appears  that,  to  a  good  approximation,  membrane 
and  curvature  effects  on  one  hand  and  flexural,  rotatory  inertia  and 
shear  deformation  effects,  on  the  other  hand,  are  mutually  exclusive 
in  two  ranges  of  wave  lengths,  separated  by  a  "transition"  wave 
length.  Thus,  in  the  full  range  of  wave  lengths,  the  associated  lowest 
phase  velocities  may  be  determined  on  the  basis  of  the  membrane  shell 
theory  (for  wave  lengths  larger  than  the  transition  wave  length)  and 
on  the  basis  of  Mindlin’s  plate  theory  (for  wave  lengths  smaller  than 
the  transition  wave  length)* 


INTRODUCTION 


Due  to  a  rapidly  increasing  technological  importance,  considerable 
attention  has  been  given  during  recent  years  to  the  analysis  of  shells  and 
the  literature  on  shell  theory  grows  at  an  accelerated  rate.  In  the  course 
of  reconsideration  and  refinement  of  the  existing  theories,  in  particular 
that  of  Love  (l)^,  various  objections  have  been  raised  as  to  the  method 
of  deriving  a  shell  theory  and  also  as  to  the  resulting  equations.  It  was 
emphasized  by  Vlasov  (2),  for  example,  that  certain  shell  theories  are  in 
contradiction  to  the  principle  of  conservation  of  energy  and,  also,  that  the 
approximations,  introduced  in  the  course  of  the  development,  are  sometimes 
dealt  with  subsequently  in  a  non-consistent  manner. 

Systematic  procedures  for  deducing  shell  theories  of  various  degree 
of  accuracy,  starting  from  three-dimensional  equations,  were  discussed  by  a 
variety  of  authors,  for  example  by  Hildebrand,  Reissner  and  Thomas  (3)  and 
by  Kennard  ( U,  5),  the  latter  basing  his  development  on  the  work  of  Epstein 
(6).  Several  such  theories  have  been  compared  with  one  another  recently  by 
Naghdi  and  Berry  (7)  in  a  qualitative  manner,  but  the  authors,  in  their  own 
judgement,  did  not  succeed  in  arriving  at  definite  conclusions  concerning 
the  relative  merits  of  the  different  approximations  used  by  various  authors 
in  the  classical  theory  of  shells. 

Thus,  the  state  of  knowledge,  briefly  described  above,  definitely 
calls  for  the  establishment  of  criteria  for  the  range  of  applicability  of 
the  various  shell  theories,  Solutions  obtained  by  different  shell  theories 

^Numbers  in  parentheses  refer  to  the  Bibliography  at  the  end  of  the  paper. 


should  be  compared  not  only  on  a  relative  basis,  but  also  with  solutions 
from  the  three-dimensional  linear  theory  of  elasticity. 

Reviewing  recently  a  controversial  paper  by  Hwang  (8)  on  the  subject 
of  shell  theories,  where  the  correctness  of  various  approaches  was  questioned, 
Truesdell  (9)  appraised  the  situation  by  stating  that  in  his  opinion  "all 
work  on  this  subject  is  purely  formal,  and  the  various  results  obtained  by 
different  perturbation  processes  cannot  be  shown  to  be  right  or  wrong  by  the 
a  priori  arguments  always  employed.  What  is  lacking  is  a  mathematical  theorem 
making  precise  the  status  of  solutions  of  any  given  set  of  proposed  equations 
with  respect  to  corresponding  solutions  of  the  three-dimensional  theory." 

The  present  paper  is  intended  to  be  a  first  contribution  in  this 
direction.  In  it,  no  precise  criteria  or  theorems  are  established!  however 
certain  quantitative  though  limited  information  is  deduced.  Comparing,  for 
the  first  time,  solutions  of  the  three-dimensional  theory  with  solutions  of 
shell  theories  it  was  natural  to  start  out  by  considering  a  simple  a  problem 
as  possible.  As  far  as  the  geometry  of  the  shell  is  concerned,  a  circular 
cylindrical  shell  was  selected,  and  since  the  authors  were  interested  in  the 
dynamic  behavior  of  the  shell,  free,  axially  symmetric  vibrations  (or  waves) 
were  studied.  Furthermore,  it  was  decided  to  test  merely  several  field 
equations,  such  that  the  shell  was  considered  as  being  infinite.  The  import¬ 
ant  and  typical  question  of  the  edge  effect  in  shell  theory  was  thus  not 
raised  at  present  at  all. 

An  analysis  of  traveling  waves  in  a  circular  cylinder  on  the  basis  of 
the  three-dimensional  theory  was  carried  out  by  Fay  (10),  but  his  results  were 
not  used,  because  it  was  noticed  that  an  error  has  slipped  into  his  development. 


Thus,  the  derivation  of  the  frequency  equation  is  presented  briefly  in  the 
first  section  of  this  paper.  A  more  complete  discussion  of  this  three- 
dimensional  theory,  including  the  transition  to  the  two  limiting  cases  of 
a  flat  plate  (Rayleigh-Lamb  solution  (11,  12))  and  of  a  solid  rod  (Poch- 
hammer  solution  (13)) ,  as  well  as  the  question  of  parasitic  modes  found  by 
Fay  (10),  will  be  taken  up  by  the  present  authors  in  a  separate  paper. 

Free  motions  in  a  cylinder  using  shell  theories  were  analyzed  by 
Rayleigh  (llj),  and  more  recently  by  Baron  and  Bleich  (1$)  and  by  Junger 
and  Rosato  (16),  but  since  only  membrane  effects  or,  in  addition,  bending 
effects  were  taken  into  account,  the  results  cannot  be  expected  to  be 
valid  for  the  full  frequency  range,  that  is  for  very  short  wave  lengths. 

Being  guided  by  the  work  of  Mindlin  (17)  in  the  theory  of  plates, 
which  includes  in  addition  to  the  classical  terms  (flexural  stiffness  and 
transverse  inertia),  terms  accounting  for  transverse  shear  deformation 
and  rotatory  inertia,  an  analogous  shell  theory  was  deduced  for  axially 
symmetric  motions  in  a  cylinder.  Thus,  in  addition  to  the  terms  bf  the 
classical  shell  theory,  (compressional  and  flexural  stiffness,  transverse 
and  longitudinal  inertia),  the  equations  derived  presently  include  the 
effect,  of  transverse  shear  deformation,  introduced  previously  into  more 
general  shell  equations  of  equilibrium  by  Hildebrand,  Reissner  and  Thomas 
(3),  and  the  effect  of  rotatory  inertia. 

The  effect  of  these  various  terms  on  the  velocity  of  traveling 
waves  in  a  cylinder  is  analyzed  and  in  each  case  the  transition  to  a 
flat  plate  is  performed.  The  "shear  deflection  coefficient,"  appearing 
in  the  relation  between  average  transverse  shear  stress  and  strain  was 
determined  in  two  ways  analogous  to  Mindlin  (17). 


Computations  of  the  lowest  wave  velocity  as  a  function  of  the  wave 
length  were  made  for  a  shell  of  thickness  to  radius  ratio  of  l/30.  It  was 
found  that  excellent  agreement  between  the  approximate  and  the  exact  solu¬ 
tion  may  be  reached  for  the  full  frequency  range,  by  including  in  a  shell 
theory  all  the  effects  mentioned  above. 

Two  fairly  distinct  types  of  shell  behavior  were  established.  For 
wave  lengths  larger  than  a  certain  "transition"  (relative)  length,  the  com- 
pressional  stiffness  is  predominant  and  the  membrane  theory  gives  good 
results.  For  wave  lengths  smaller  than  this  transition  length,  the  shell 
behaves  essentially  like  a  plate  and,  for  good  results,  calls  for  inclusion 
of  transverse  shear  deformation  and  rotatory  inertia  effects.  Thus,  the 
use  of  bending  theory  of  shells  appears  to  be  of  little  value.  In  the  full 
frequency  range,  vibrations  of  a  shell  may  be  described  with  good  accuracy 
by  the  membrane  theory  and  by  Mindlin's  plate  theory.  It  may  be  expected, 
that  this  statement  will  not  be  valid  for  non-axially  symmetric  motions, 
which  will  be  studied  by  the  present  authors  in  a  forthcoming  paper. 

Since  interest  concentrated  at  present  on  free  vibration  solutions, 
no  attention  was  given  to  terms  arising  from  external  loading  and  no  energy 
functions  were  constructed,  which  would  permit  the  establishement  of 
appropriate  initial  and  boundary  conditions,  sufficient  to  assure  a  unique 
solution.  Furthermore,  attention  was  restricted  to  the  lower  mode  only  and 
no  mode  shapes  were  investigated.  It  is  intended,  however,  to  study  these 
questions,  when  dealing  with  a  more  general  theory  of  cylindrical  shells, 
valid  for  non  axially-syiranetric  motions. 
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SQLUTION  BT  THE  THREE-DIMENSIONAL  THEORY 

The  circular  cylindrical  shell  of  inner  radius  a  and  outer  radius  b, 
is  referred  to  cylindrical  coordinates  r  Br%  Xf  the  x-axis  coinciding  with 
the  axis  of  the  shell.  Designating  by  Lf  U  U,the  components  of  the  dis- 
placement  in  this  coordinate  system*  attention  will  be  given  to  axially 
symmetric  motions  only*  that  is  to 

“«■  * 0  ■  h  * 0  111 

The  linear  elasticity  equations  of  motion,  for  a  homogeneous*  isotropic  shell 
material  obeying  Hooke's  law*  may  then  be  given  in  the  form 

^  +  [2] 

-  2/1  -puK 

where  the  dilation  A  is 

A  =  i  Xrur)  t  [3] 

r  $r  ix 

and  the  rotation 


Aj^U^re  Lame's  constants  and  p  is  the  mass  density.  Dots  indicate 
differentiation  with  respect  to  time*  t  . 

Since  the  shell  is  assumed  to  be  free  from  surface  tractions*  the 
boundary  conditions*  to  be  imposed  on  the  cylindrical  surface  of  the  shell, 
are 

For  r  =  a)  b  :  <£r  =  =  o 


[5] 


In  solving  the  equations  of  motion  [2],  it  is  found  more  convenient 
to  consider  A  and  as  the  dependent  variables,  rather  than  C4randUx» 
The  equations  [2]  take  then  the  form 


it  + 1  w  f  ^ 

IP  r  ir  ax* 


=  P  >*4 


it 


A+ijt  it 

—t  +  -  3  -  Q  +  i\  =  P  i'tOf 

it  t  p»  jx1 

Assuming  the  motion  to  be  harmonic 

i(  tot-xx) 

%  = 


[6] 


[7] 


where  U)  is  the  circular  frequency  and  the  phase  function,  related  to 


the  phase  velocity  C  by 

eX  =  ^/C 

[8] 

the  equations  of  motion  [6]  become 

/a 

TP 

t ■  l  i&  +  &XA  r  0 

t  IF  ' 

[9] 

iP 

+.  1  ^  H  r  =  o 

r  il"  "P 

where 

r  c;  ~  c') 

y1  =  "Yp"  £>) 

vV  V" 

[10] 

and  Cc  ^  Cs  are  the  phase  velocities  of  congressional  and  shear  waves, 


respectively 
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Equations  [9]  are  of  the  Bessel  type,  and  A  and  ^must  be  taken  proper- 
tional  to  the  Bessel  functions  7^))  7(\f)  )  Y  (ft) 

respectively.  The  displacements  Ur>Ux  themselves,  in  view  of  the  relations 
^  [li]>  become  therefore 

>re  A,  B,  C  and  are  arbitrary  constants  of  integration. 

The  boundary  conditions  [5]  in  terms  of  the  displacements  will  be 

For  r  =  Oj  b  : 

A  A  +  i/*-  _  Q 

/<■(!£+  &*)  =  o  tl3] 

'  ?X  >|r*  / 

These  four  boundary  conditions  yield  four  homogeneous,  algebraic  equations 

in  the  constants  A,  B,  C  and  D.  For  a  non-trivial  solution,  the  vanishing 

teir  determinant  is  required,  which  results  in  the  characteristic 
equation 

HK)  *<fh'fnb) 

-  (l+B)/xt>b][aKJfl  t 

-[(i+i)  /Fnbja  K„  (()  +  6  K ,  (r)]  K  (ft 

-  O 


[11*] 
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where 

« mW  =  IMXM  - 


It  may  be  of  interest  to  note  that  the  last  three  terms  of  the  characteristic 
equation  [lU]  appear  to  have  been  omitted  by  Fay(L0), 

For  computational  purposes,  the  equation  [lU]  is  written  conveniently 

in  two  different  forms,  depending  upon  —  ^  I  <,  In  the  range  £.  <  | 

cs  Cg 

the  arguments  of  the  Bessel  functions  become  imaginary  and  equation  [lU] 
takes  the  form 

U")  =  [ tMJpMJ?)  - 

-  [(i+S)/tofc  »J}>  -  k  MJfllMJi)  1161 

i-[(n-eJ/Ffafc][a  M Jr)  - 

=  0 
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where 

^  r  -/S*  [17] 

Y*  =  -  rl 

F  =  sV/^'F 


I  )  and  /{  are  the  Modified  Bessel  Functions  of  the  first  and  second 
kind,  in  the  notation  of  Watson  (18). 


In  the  range  £.  <  £.  ,  fi  is  imaginary  and  Y  is  real. 

9 $  * 

Equation  [lli]  may  then  be  written  in  the  form 


=[W  -MJ)K, .(*)+(*/, r/irok) 

-[(!+e)/nb][aM,.(s)  - 1  MJfijK/i) 

-  [(/f  eVfrok  ][«  *,.M  +  !>  KJri]M {[}) 

'Ib+V/rt'ckWfK  W 


where 


[19] 


Solving  equations  [16]  and  [18],  respectively,  the  phase  velocity  $  :  £ 
r  , 

versus  o  r  h  (with  VrO'3*  !l  r—  )  is  plotted  as  curve  (T)  in  Fig.  1. 

Z.  *  30  ^ 

For  sufficiently  thin  shells,  the  arguments  of  the  Bessel  functions 

appearing  in  the  characteristic  equation  [16]  become  very  large  and  a  simpli¬ 
fication  can  be  made  by  expanding  these  functions  asymptotically.  In  this 
case  the  frequency  equation  appears  in  the  much  simpler  form: 


2  £  cosh  cosh  ?h  —  /  j 
-  (F+-  zr)sinhfih  sinhyh 

+  [(!  +  *)* /Fr'cib]  Sihhjih  sihbfh  [20] 

+-  [(H-B)h  /fab]  Sihh  Vhcoshfh 

~  [(it  B)h  /f  yab]  s]r'hfb  coshir  h  -o 

where  h  is  the  thickness  of  the  shell, 


Limiting  Velocities 


For  very  short  wave  lengths,  h*-**,  thusyJ-*00*  and  for  very 

large  arguments  sinh  z  ~cosh  z.  Hence  the  frequency  equation  [20]  takes 
the  form 

^  -  (F+j)  =  O 


that  is 


F  =  / 


or 


J(nl-sl)(l-s)  =  ri(±S*- if  o<S<  I 


[21] 


where 


s  =  £ 


5 

n1  =  £  s  z(i-»)/(i-2»)  [22] 


Expression  [21]  is  the  equation  for  the  velocity  of  Rayleigh  surface  waves. 
Thus,  very  short  waves  are  propagated  in  the  shell  with  the  same  velocity  a 
in  a  plate,  the  latter  case  having  been  investigated  by  Lamb  (12). 


{I 
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For  very  long  wave  lengths,  h  «  O  and  hence,  in  the  limit, 

O  •  It  may  be  seen  that  in  this  case  —  >  |  .  The  arguments 
of  the  Bessel  functions  become  then  very  large  and  on  expressing  these 
functions  again  in  asymptotic  form,  the  frequency  equation  [llj]  reduces  in 
the  limit  to 


<C 

)-iV  =  0  £  >  / 

cs 

ft 

or  C 

cs 

=  2  (l+v) 

[23] 

ft 

that  is  q 

=  £/? 

which  is  the  "bar"  velocity,  as  it  is  well  known  from  elementary  bar  theory 
and  as  it  was  determined  by  Pochhammer  (13)  for  the  analogous  limiting  case 
in  his  study  of  motions  of  solid  cylinders. 

Thickness-Shear  Motion 


•>  h  -  ° 


If  the  displacement  components,  in  addition  to  assumptions  [1],  are 
restricted  still  further  by  requiring 

I2h) 

the  possible  motion  is  designated  as  thickness-shear  motion  (17).  The 
equations  of  motion  [2]  simplify  then  to 

c*[  +  i  K 1  _  $ 

s  L  r  ip  J  "  * 


3i-‘  r  ir 

Assuming  the  motion  to  be  periodic 

=  F(ir)e 


[25] 


ciot 


[26] 


equation  [25]  reduces  to 
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1  dF  +  «£  F  _  o 
1 'dr  c? 


[27] 


which  is  a  Bessel  equation  of  aero  order.  Imposing  boundary  conditions 
analogous  to  [13]  (traction-free  surfaces),  the  resulting  frequency 
equation  may  be  put  into  the  form 


X&*)  ’  yw 

4  =  0/t 

I  =  ui,lct 

Thus,  for  a  given  value  of  -A  ,  I  can  be  determined  and  hence  & . 


[28] 


[29] 


For  thin  shells,  the  arguments  of  Bessel  functions  in 
[28]  become  large  and  the  frequency  may  be  computed  from  th  simplified 

equation 

ton  =  °  [30] 


which  was  obtained  after  representing  the  Bessel  functions  by  asymptotic 
expansions.  The  lowest  frequency  from  [30]  is 

CO  -  T"  Si  [31] 

h 

which  is  identical  to  the  one  obtained  in  studying  the  corresponding 
motion  of  a  flat  plate,  see  (12)  and  also  (17). 
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APPROHMATE  SHELL  THEORIES 


In  the  theory  of  shells ,  it  is  customary  to  employ  t  coordinate  system, 
fomed  by  two  axes  being  in  the  middle  surface  of  the  shell  and  the  third 
axis  being  normal  to  the  middle  surface.  Thus,  the  coordinates  will  now  be 
taken  as  0j  2 )  X  being  still  in  the  axial  direction  of  the  cylindrical 
shell  and  the  new  incoordinate  being 

2  -  r-R  [32] 

where  R  -  Cl  +  b  is  the  mean  radius  of  the  shell. 

2 


Stress-Equations  of  Motion 


Por  the  case  of  axial  symmetry,  the  stress  equations  of  motion  may 
then  be  written  as 


+ 

IF 

■h  Jk  =f>iiK 

[33] 

+■  ^  Oxl 

+ 

fu 

l 

II 

“D 

Ik 

— - — ;  )  * 

R  +  f 

It  is  now  assumed,  that  the  displacements  Ux  and  of  the  three- 
dimensional  theory  may  be  approximated  by  displacements  U  and  U  ,  whose 

*  c 

dependence  through  the  thickness  h  of  the  shell  shall  be  specified  in  the 


form 

Ux  -  u(X>  t)  t  2  %  (x,  t) 
ug  =  ur(x,t) 

The  shell  stress-equations  of  motion  are  obtained  by  substituting  in 


[3U] 


equations  [33]  U„ ,  U.  by  .  U_  ,  respectively,  multiplying  the 

*  i  e  *  /  t 

first  and  the  second  by  (R  +  3)  and  the  second  of  equations  [33] 

and  integrating  through  the  shell  thickness  h  ,  from  3  =  - 1  to  ?=  h. 

*  2 
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The  justification  for  the  multiplication  mentioned  above  stems  from  energy 
considerations  which  shall  be  discussed  by  the  authors  elsehvere. 


The  resulting  equations  are: 

mJ  -  0„  *  piYvi  +  ?«) 


V= 


where  use  has  been  made  of  the  fact  that  the  shell  is  free  from  surface 

tractions,  i.e.  r  <T  =  O  on  £  -  £  H  .  Primes  indicate  differ- 

**  i»i  J 

entiation  with  respect  to  x. 

The  stress  resultants  (or  shell  stresses)  are  defined  by  the 


following  integrals: 


* 

-  J  ‘  * 

-.S'  r"0,Vldl 


In  case  of  equilibrium,  i.e.  if  all  the  time  derivatives  vanish, 
the  equations  [35]  represent  a  special  case  of  a  more  general  shell  theory, 
deduced  by  Hildebrand,  Reissner  and  Thomas  (3).  If  the  radius  (of  curvature) 

R  tends  to  infinity,  that  is,  the  shell  approaches  a  flat  plate,  the  first 
two  equations  of  the  set  [35]  reduce  to  those  derived  by  Mindlin  (17 ),  which 
govern  flexural  motions  of  plates,  including  the  influence  of  rotatory  inertia 
and  shear.  The  third  equation  represents  longitudinal  motions  of  the  plate, 
which  are  not  coupled  to  flexural  motions.  The  terms  with  ph  represent  the 
influence  of  rotatory  inertia.  It  is  seen,  that  there  are  altogether  three 
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such  tens.  The  ter*  misrepresents  transverse  inertia  and  the  tens  phu 
is  longitudinal  inertia. 


Stress-Displacement  Relations 


The  relation  between.  the  shell  stresses  defined  by  the  set  [36] 
and  shell  displacements  introduced  by  [3U]>  will  be  derived,  starting  from 
the  stress-strain  relations  of  an  isotropic  three-dimensional  solid  obeying 
Hooke's  lav.  If  the  relation  for  (T,  is  solved  for  0  and  then  substituted 
into  the  relabiomfor  the  other  two  normal  stress  components,  the  result  is 


where  £  is  Young's  modulus  and  V  Poisson's  ratio.  This  procedure  is 
justified  by  observing  that  would  vanish  on  the  basis  of  the  form  given 
to  the  transverse  displacement  On  the  other  hand,  unrestrained  effect 

of  Poisson's  ratio  will  take  place  in  the  actual  shell  in  the  z-direction. 
These  two  mutually  exclusive  requirements  may  be  reconciled  by  eliminating 
€^2  from  the  stress-strain  relations  in  the  manner  described  above. 

In  addition  to  the  normal  stresses,  the  shear  stress  (FJ^is  needed, 
which  is  related  to  the  shear  deformation  Y  by 

^ X2  =  ^  ^X2  ^8] 


where  G  is  the  shear  modulus.  With  the  aid  of  the  strain  displacement 


relations 


= 

ix 

=  u*  l(R+i) 

_  H:  +  iWj 

"  n  u 


[39] 
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th«  approximate  strains  are  evaluated  using  assumptions  [3b]  and  substituted 
into  the  stress-strain  relations  [37],  [38].  Then,  integrations  are  carried 
out  as  specified  by  the  forms  [36].  The  results  are  modified  in  three 
respects t 

1.  The  integrals  containing  OT^are  omitted,  just  as  in  the  classical 
beam  or  plate  theories. 

2.  the  Integra:  containing  Yx2  is  multiplied  by  a  constant 
(the  shear  coefficient),  to  be  determined  later. 

3.  logarithmic  terns  occuring  in  the  integration  are  in 

series  of  ascending  powers  of  h  and  terms  higher  than  the  cubic  one  are 

K 

neglected. 


The  resulting  shell  stress-displacement  relations  are 


equations  of  motion,  displacement  equations  of  motion  are  obtained,  which 
may  be  put  in  the  form 
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-[ K*S k  jL  k  +[£  *  £  =  o 

^f,]- 

s  *  [1*2] 

I  is  the  moment  of  inertia  h.  ,and  £  is  the  plate  compressional  modulus 

f, «  « 

This  form  Is  convenient  in  the  discussion  of  the  physical  significance 
of  the  terms  influencing  the  motion n  There  are  three  stiffnesses  involved,  D 
representing  the  influence  of  flexural,  £p  of  compressional  and  G  of  trans¬ 
verse  shear  stiffness,  respectively.  In  addition,  there  are  three  inertia 
terms  involved,  pi  representing  the  influence  of  rotatory  inertia,  ph  in  the 
second  equation  representing  transverse  inertia  and  ph  in  the  third  equation 
representing  longitudinal  inertia.  The  rotation  is  coupled  to  the  trans¬ 
verse  displacement  UT  through  the  transverse  shear  stiffness  and  to  the  long¬ 
itudinal  displacement  U  through  flexural  stiffness  and  rotatory  inertia. 

The  coupling  between  displacements  L4  and  UT  is,  essentially,  both  a 
Poisson's  ratio,  and  a  curvature  effect. 


For  a  thin  shell,  terms  with  V/R  and  I fft  may  be  omitted,  because 
D  and  I  are  proportional  to  the  cube  of  the  shell  thickness  h  •  As  is  seen 
from  equations  [1*2],  the  suppression  of  influence  of  flexural  stiffness  and 
rotatory  inertia  implies  the  assumption  of  a  thin  shell,  but  not  vice  versa, 
because  in  the  first  bracket  of  the  first  equation  [1*2],  D-and  I- terms  are 
unaffected  by  the  curvature  .  For  a  thin  shell  we  have  then 


-18- 


C®n« -*’5h ~  £ ’l‘5,,|,]“r  = 

[ -  K‘S  H  h  +[%-*!•&  Th  |>  >.  Ju  . 

[%Utr  + &  55‘- '  P'1  ftJ‘ «  = 

and  there  is  no  direct  coupling  between  y>  and  d  .  Doing  one  step  further^ 
the  shell  can  be  Bade  to  approach  a  flat  plate,  that  is  R-+  oo  ,  The  first 


two  equations  reduce  then  to  Hindi in's  piste  theory  (1?)  and  the  last,  not 

coupled  to  the  first  two,  ie  the  classical  equation  for  compressional  motions 
in  a  plate. 


Wishing  to  suppress,  in  the  more  complete  equations  [k2],  the  influence 
of  transverse  shear,  we  may  not  set  the  terms  with  G  simply  equal  to  zero  in 
)se  equations,  because  these  terms  represent  the  transverse  shear  force, 

while  we  wish  to  neglect  transverse  shear  deformation,  which  is  effected  by 
setting 

%  =  m 


There  be  only  two  displacement  components  left  and  the  terms  with 

G  and  may  be  eliminated  with  the  aid  of  the  above  relation  [hk],  The 
result  is 


f 
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a  «.  „  .to.  *  ^  ^  ^ 

shells,  see  for  example  Flugge  (20)  or  Vlasov  (2). 

For  thin  shells,  J>/R=  I/R  s0f  they  peduce  ^ 

f +  5  ^4  -pi )*  7  [£v  i  7 

K1  E*  SxTt'--'^  Haifjueo 

*  R  >*  tb6] 

[  —  ]  ur  f  f  f  >2  -I  .i  1 

*  >*J  ^Tx'  ~  ~° 

It  the  rotatory  inertia  ten,  is  neglected  and  the  transition  to  a 
flat  plate  is  performed  ^  ,  the  two  motions  in  ^  and  U  becc*. 

*g8in  'C0Upled-  Th6  "otion  of  Ur  will  be  governed  by  the  classical  bending 

o  of  plates  and  the  motion  of  U  will  be«ain  that  of  classical  com- 

pressional  theory. 

If,  instead,  the  rotatory  inertia  term  and  the  flexural  stiffness 
term  are  neglected  in  equations  (1*6),  the  result  is 

4  ,o 

Which  represents  the  familiar  membrane  theory  of  shells. 


WAVE  SOLUTIONS 

We  consider  now  the  propagation  of  free  harmonic  waves  and  seek 
solutions  of  shell  equations  in  the  form 


=  fQi&t-**) 
ur  (x,  t)  =  We'^'^ 


I 
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Design&ting  by  I  the  wave  length,  the  natural  frequency  (O  and  the  wave 
number  o<  ,  are  related  to  the  phase  velocity  C  by 

CO  =  2fC  •  oC 

L  ' 

Introducing,  further,  the  notation 

m  =  h/R 

S  =  C/Cs 

f  =  h/L 

N  ~  1 

and  substituting  forms  [!|8]  into  the  most  complete  shell  equations  [1*2]  of 
the  present  investigation,  the  characteristic  equation  is  calculated  to  be 

i  !2i  ) 

3  L  12  J  12 

-4k*sV*  [51J 

-  *  jjV/v] 

3  J 

—  N  V*  j  =  O 

IP  J 
2 

There  are  thus  three  roots  for  $  ,  but  attention  will  be  given  only 
to  the  smallest  one.  In  the  whole  range  0<  the  values  obtained 

differed  too  little  from  those  obtained  by  the  solution  of  three-dimensional 
equations,  to  be  discerned  on  the  scale  of  Fig.  1.  Thus,  curve  Qin  Fig.  1 
represents  also  a  solution  of  equation  [51].  Numerically  largest  discrepancies 
between  the  two  solutions  were  obtained  in  the  neighborhood  of  the  minimum 
value,  which  may  be  due  to  a  coupling  effect  with  purely  radial  vibrations. 

The  limiting  velocities  may  be  evaluated  directly  from  equation  [5l]. 

For  very  long  waves,  the  limiting  velocity  is 

Sl  =  [iO*-*)  +  Hl'/ifr-v)]  fO^—1 ) 


2W 

L 


[1*9] 


[50] 
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If  ■  is  small  as  compared  to  unity* 

S  =  2(j+») 

which  coincides  with  the  exact  solution*  given  by  equation  [23]. 

For  very  short  waves*  f-^oo  9  equation  [51]  reduces  to 

(2N-S,)l(K- s‘;  =0  m 

whatever  the  shell  parameter  m. 

The  double  root  $  r  2N belongs  to  the  two  higher  modes  and  will  not 


be  discussed  presently.  The  root  S*s  K* belongs  to  the  lower  mode  and  the 
coefficient  X  may  be  now  determined,  precisely  as  in  Mindlin's  plate  theory 
(17),  that  is  by  matching  the  velocity  obtained  above*  with  the  velocity 
obtained  from  a  solution  of  the  three-dimensional  theory,  given  by  equation 
[21].  An  appropriate  value  of  K  is  therefore  the  lowest  root  of 


J  0~  *')(*'- K*)  =  0<K<\ 

Thus*  K  depends  upon  Poisson’s  ratio  and  ranges  from  0*76  for 

VzOto  X  *  =  0*9/  for  V  r  0-5.  For  ^  =  0«3,  X*  =  0'M. 

Letting  m-*0,  the  shell  approaches  a  flat  plate  and  equation  [51] 
reduces  in  part,  to  the  one  obtained  by  Mindlin  (17), 


Y  V- s'X?  -  &  =  i 


which  gives  for  long  wave  lengths,  S -f  O 

si  = 

3 
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}  pi 

that  is  $  =  O  for  0  =  O ,  and  yields  thus  a  different  limiting  value  than 
for  the  shell*  For  short  wave  lengths,  ,  the  shell  equation  [51]  end 

the  plate  equation  [51*] ,  reduce  to  the  same  expression,  equation  [52]. 


Neglecting  Rotatory  Inertia 


Omitting  now  the  terms  with  pi  in  equations  of  motion  [1*2],  the 
wave  solutions  [1*8]  lead  to  the  characteristic  equation, 

(in- s\)[  l  w  -  IN*1 MV  -  ♦ffVsV1 

+  4  nW  mY,+—  )  + 

*3  /2  [55] 

+  [  ^A/VW^S-x1)  -  itfVfVr) *(»+£*) 

-  ?  w'vxVV’m1  -  4  -  r  O 

3  3  J 


The  solution  of  equation  [55]  is  plotted  as  curve  @  in  Fig.  1.  For  long 
waves,  *"-> O  ,  and  far  short  waves,  S  »  X  .  For  a 

flat  plate,  fn  -  0  ,  and  equation  [55]  reduces,  in  part,  again  to  the 
one  discussed  by  Mindlin  (17),  in  the  corresponding  case 


s’  =  ii!£' f  / 1- f‘J 

3/i-y)  L 


Neglecting  Shear 

Using  equations  of  motion  [1*5],  the' membrane  effect,  flexural  stiffness 
and  rotatory  inertia  are  included*  The  motion  will  consist  of  just  two  modes 
and  the  characteristic  equation  is  obtained 
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(2N-s‘)[-S  rV*  +l*'NSV  +  *  tV/’ 

3  f  3 

-  ±  rW$*f*  +  4t r*$  V**+  i 


-J/W/*V)7 

12  yj 


157] 


+  ^y/V'fr)1  r  1  TT^A/;1  _  T*SV*  +  iV/v7  _ 

3  y  J 


The  solution  of  equation  [57]  is  plotted  as  curve  @  in  Pig.  1.  For 
c  wave  lengths,  J""->  O  ,  the  lower  velocity  is 

S*  =  [s(i+v)+  ™'/i  ()-»)]  /(li-HL) 

For  short  wave  lengths,  f -»  oo  ,  the  lower  waye  yeloclty  ^  be  giyen 
by 

sl  = 

i-y 

Which  is  identical  to  the  one  obtained  by  Mindlin  (17)  and  i.  much  too 
large  as  compared  to  the  Rayleigh  wave  Telocity.  If  m  =  O  ,  Mindlin'a 
(17)  corresponding  result  for  the  plate  is  obtained 

s*  =  jiLrr,+zlrr' 

L  3  J  [58] 


Bending  Theory 

The  bending  theory  of  shells,  equations  [1,5],  with  1=0  leads  to 
the  characteristic  equation 

(*N'S -Ia/ttVM  m 

+  4A/‘w‘  +  ifj1)1  =  c 

3 
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The  solution  of  equation  [59]  is  plotted  as  curve  0  i»  fig.  i.  Kor  long 
wave  lengths,  ,  the  lower  Uniting  velocity  is 


For  short  wave  lengths,  /-*  ee  , 

S  -*  S.  -  lr> 

With  m=0,  equation  [59]  siaplifies  to 

sl  =  jilr1 

3  (i-») 

which  is  a  result  of  the  classical  plate  theory. 


Membrane  Theory 

The  classical  membrane  theory  of  shells,  equations  [J»7]  gives  rise 
to  the  characteristic  equation 

2TT  $  $  —  b^TTSf  —  NS^in*  +  s;  O  ^ 

which  still  contains  two  modes e 

The  solution  of  equation  [61]  is  plotted  as  curve  0  in  Fig.  l. 
For  long  waves,  f-*0  ,  the  lower  mode  has  a  velocity  S  *  =  tQ+v) 

For  short  waves,  f-*co  ,  the  lower  mode  has  a  vanishing  velocity.  For  a 
flat  plate,  rtl  =  O  ,  equation  [61]  reduces  to,  for  any  wave  length, 

Sl  =  3  N 

or 

c8r  ef p (,-„>) 

which  is  the  compressional  wave  velocity  in  a  plate.  This  is  the  only  mode 
propagated  with  a  finite  velocity.  It  is,  however,  the  mode  which  for  mjo 
is  associated  with  a  higher  wave  velocity. 


Thickness  -  Shear  Motion 

As  an  alternate  way  of  determining  the  constant  X2  in  plate  theory, 
Mindlin  (17)  considered  the  thickness-shear  vibration  and  determined  K* 
by  matching  frequencies  obtained  by  exact  and  approximate  theories. 
Precisely  the  same  can  be  done  for  a  cylindrical  shell. 

Setting  in  the  present  theory  ;  (c>  t 

u  =  ur  r  o  ^  =  Ae 

in  equations  of  motion  [1*2],  the  frequency  00  is  obtained  as 


which  is  independent  of  the  radius  of  curvature  R  • 

Equating  the  right-hand  sides  of  equations  [31]  and  [63]  results  in 

KX  =  5  [61*] 

12 

The  difference  between  K  given  by  equations  [53]  and  [6U]  is  not  large 
and  vanishes  if  Poisson's  ratio  V  ■  0.176. 

Discussion  of  Results 

The  essential  results  of  the  present  study  are  visible  from  Pig.  1. 
The  important  conclusion  which  can  be  drawn  is  to  the  effect  that  it  is 
possible  to  deduce  a  relatively  simple  shell  theory,  which  predicts  for  the 
full  range  of  wave  lengths  practically  the  same  phase  velocities  of  prop¬ 
agated  waves  (or  frequencies  of  free  vibration),  as  the  three-dimensional 
theory. 

Furthermore,  it  is  seen,  that  for  large  wave  lengths,  £ <  0.035*  the 
elementary  membrane  shell  theory  yields,  as  expected,  good  results.  From 


this  value  on,  that  is  if  the  wave  length  L  is  smaller  than  the  mean 

radius  of  the  shell  R,  (since  ft =?  £l  E  and  h  was  taken  as  y-  )  the 

R  L  J  3° 

shell  bending  theory  must  be  used,  which  predicts  a  minimum  of  the  phase 
velocity.  An  approximate  value  for  this  minimum  may  be  obtained  by  ob¬ 
serving  that  since  bending  is  important,  in  this  region,  the  effect  of 


axial  inertia  may  be  neglected.  The  minimum  velocity  is  thenfound  to  be 

s‘  = 

"  IfZW 

and  the  corresponding  value  of  is 

£*  =  "i nEEi.  ^ 

T  ' 

which  gives 


-  no 


The  expression  for  the  minimum  velocity  suggested  by  Junger  and  Rosato  (16), 
differs  markedly  from  the  one  given  by  equation  [65].  For  wave  lengths 
larger  than  those  given  by  equation  [66],  that  is  if  1  is  smaller,  than, 
say  twice  the  radius  R,  the  curvature  terms  are  practically  of  no  importance 
and  the  shell  behaves  like  a  flat  plate,  but  rotatory  inertia  and,  much  more 
so,  transverse  shear  deformation  effects,  become  very  significant. 


In  summarizing,  it  may  be  stated  that  a  thin,  or  medium  thick,  shell, 
behaves  radically  different  in  propagating  a  harmonic  wave,  depending  upon 
whether  the  wave  length  is  smaller  or  larger  than  the  "transition"  length 
given  by  equation  [67]  0  For  wave  lengths  larger  than  the  transition  length, 
the  shell  may  be  well  described  by  the  membrane  theory,  because  all  bending 
effects  are  negligible  and  curvature  terms  are  predominant.  For  wave  lengths 
smaller  than  the  transition  length,  the  shell  behaves  like  a  flat  plate,  since 
curvature  terms  become  negligible.  Thus  it  appears  that  actually  in  the 
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problem  on  hand,  no  advanced  shell  theory  is  necessary.  The  full  range 
may  be  described  by  the  membrane  shell  theory  (for  large  wave  lengths)  and 
by  the  Mindlin  plate  theory  (for  short  wave  lengths). 

The  existence  of  a  minimum  may  be  made  plausible  physically  by 

/ 

observing  that  as  the  wave  length  becomes  smaller,  the  contribution  of 
radial  displacements  to  the  motion  becomes  larger  and  therefore  the  coupling 
with  purely  radial  vibrations  becomes  stronger,  exhibiting  in  the  minimum  a 
certain  resonance  effect.  After  a  further  decrease  of  wave  length,  the 
motion  becomes  essentially  flexural  and  the  coupling  gets  weaker.  A  certain 
improvement  of  the  shell  theory  offered  in  this  paper  may  be  achieved,  by 
matching  the  frequency  of  purely  radial  vibrations  obtained  from  present 
theory,  with  the  one  determined  from  three-dimensional  theory  (21).  This 
would  necessitate  the  appropriate  introduction  of  an  additional  constant  in 
the  stress-displacement  relations  [1*0]  of  the  present  shell  theory,  which 

i 

would  play  a  role  similar  to  the  constant  X  • 

In  conclusion,  it  should  be  emphasised  that  the  above  remarks  are 
valid  only  for  lower  phase  velocities  and  are  not  necessarily  true  for 
velocities  of  higher  modes,  for  mode  shapes,  for  stresses  in  the  shell  and 
for  group  velocities.  Neither  are  they  necessarily  valid  for  non  axially- 
symmetric  motions,  which  shall  be  investigated  by  the  authors  in  a  forth¬ 
coming  paper. 


- 
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